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Abstract 

In this paper we suppose that the cosmological constant will change when the universe expends. For a 
general consideration, the cosmological constant is assumed to be a function of scale factor and Hubble 
constant. According to the ADM formulation, to the FRW metric, the extrinsic curvature / equals -6H 2 
and spatial curvature R equals 6k /a 2 . Therefore we suppose cosmological constant is a function of extrin- 
sic curvature and spatial curvature. We investigate the cosmological evolution of FRW universe in these 
models. At last we investigate two particular models which could fit the observation data about dark en- 
ergy well. Actually a changeless cosmological constant is not essential. If when the universe expands, the 
cosmological constant changes slowly and gradually flows to a constant, the observation data about dark 
energy could also be fitted well by this kind of theory. 
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I. INTRODUCTION 



Eintein's general relativity(GR) has been considered as a fundamental theory of gravity. Even 
from an effective field theory point of view, it should describe all large scale gravitational physics, 
in particular the evolution of our universe. However, the discovery of dark matter and dark energy 
from various observation posed great challenge to the theory. To address the dark energy issues, the 
most simplest method is adding a cosmological constant which is called ACDM model. The cos- 
mological constant problem has been investigated by many papers such as [0-0] and so on. Many 
other theories also have been pro pos ed for the cosmic acceleration, for instance, Quintessence[|5|- 
0], K-essence[[8|,|9l], PhantomlflO .Hi 1 . etc. And numerous versions of modification or extension of 
GR have been proposed in the last and this century. The /(i? (4) ) theory Ill2l - ll5n is one of the modi- 
fications to GR. In this theory, the Lagrangian density / is an arbitrary function of i? (4) . The theory 
can explain the cosmic acceleration without introducing a cosmological constant, even though the 
theory would be very much restricted by the solar system test. 

In this paper, we suppose the cosmological constant will change when the universe expends. 
For the scalar factor and Hubble constant are main variable to characterize a universe, the cos- 
mological constant is supposed to be a function of them. In the ADM formulation yj], the metric 
is 

ds\ = -N 2 dt 2 + gijidx 1 - N l dt)(dj - N l dt). ( 1 ) 

the building block to construct the action is the gauge invariant quantity: the extrinsic curvature 
tensor 

Kij = ^ (gtj ~ ViNj - VjN t ) , (2) 

where a dot denotes a derivative with respect t, and the 3-dimensional spatial curvature R comes 
from gij. In terms of them, the action of GR with a cosmological constant is 



EH 



= J dtd 3 x [KijK iJ - K 2 + R - 2A] 



(3) 



where K = g' J Kjj. And we define extrinsic curvature / as 

/ = KijK ij - K 2 . (4) 
The homogeneous and isotropic universe is described by the FRW metric, 

ds 2 = -dt 2 + a(t) 2 | - dr k ^ + r 2 (d6 2 + sin 2 6># 2 ) j , (5) 

where k = 0, +1. To this FRW metric, After some trivial calculations, it is easy to get the extrinsic 
curvature / = -6H 2 and spatial curvature R = 6k/a 2 . Obviously / is proportional to the square 
of Hubble constant and R is proportional to the inverse square on scalar factor a. Therefore we 
suppose the cosmological constant A = A(7, R) as a function of extrinsic curvature / and spatial 
curvature R. 

For the calculational convenience, we define F(I,R) = I + R - 2A(I,R) and the action of this 
theory is 

S P = J dtd 3 x ^NF (/, R) (6) 

This theory could be regard as an extension of f(R^ A) ) gravity model. In the next section, we 
exhibit the basic equations by varying the action with respect to the function N,Nk,gij. Then we 
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investigate the cosmological evolution of FRW universe in these models. At last, we investigate 
two particular models which could fit the observation data about dark energy well. Actually a 
changeless cosmological constant is not essential. If the cosmological constant changes slowly 
and gradually flows to a constant, the observation data about dark energy could also be fitted well 
by this kind of theory. 

H. BASIC EQUATIONS 

Let us start from the action 



Actually the action could be more general if F is an arbitrary function of " gij,K,Kij, 
ViK jk , • • • , V h V h ■ ■ ■ V in K jk , ■ ■ ■ ,R,Rij,R ijkU ViR jk , • • • ". In this paper we would like to focus on 
the simplest case F = F (I,R). Now we consider the case that there are matters couple to the 
gravity. The action should be 



By varying the action with respect to the function N,N k ,gij, we get the Hamiltonian constraint, 
super-momentum constraint and dynamical equations, 




(7) 




(8) 



1 



[F-2FjI], 



(9) 



\6nG 
1 _ 



8ttG 



(10) 




(ID 



n lj in above equations is defined as n' j = K lj - AKg lj . And J N , J k , p lj are defined as 




(12) 




(13) 




(14) 



(15) 
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III. COSMOLOGICAL MODELS 



The homogeneous and isotropic universe is described by the FRW metric, 

dr 2 



ds 2 = -c 2 dt 2 + a(t) 2 



l-kr 2 

where k = 0, + 1 . The extrinsic curvature of the spatial slices is 



+ r 2 (d6 2 + sin 2 &/0 2 ) j , 



(16) 



_ i«(0 

c a(t) 



K = 1 3a(Q 
c a(t) 



(17) 



Then / = -6H 2 where H = a/ a which is defined as Hubble constant. The spatial slices are of 
constant-curvature 

2k 6k 

R " = sm g * R = WY <18) 

It can be shown that the momentum constraint (fTOl) is satisfied identically when J k = for the 
perfect fluid. From the Hamiltonian constraint ® and the dynamical equation (fTTI) . and "c" is 
rescaled to unity, we get two cosmological equations as 



F+ 12F (/ (-j = \6nGp m , 

ld\ 2 2k 1 
+ 4Fj - - F r— + -F = -SnGp m , 
\al a 1 2 



„ a .a 
Fj- + Fj- 

a a 



with the energy conservation equation 



p m + 3-(p m + p m ) = 0. 
a 



(19) 
(20) 

(21) 



The FRW equations with the matter and the dark energy could be rewritten as 

(Pm + PDE) , 



d\ 2 k _ SnG 

a) a 2 3 

2 k 



(22) 
(23) 



Here we have regarded the nonlinear terms in the equations (1191.(1201) as the dark energy effectively. 
Comparing ((221) to (HU), ([23]) to dH, it is easy to get 



a (a\ k „ 
2- + - + — = -SnG(p m + p DE ) . 
a \al 



SnGp DE = 3(1- 2F d ) (-) + ^-^F, 

\al a L 2 



SnGp DE = 2 (Fj - 1) - + 2Fj- + (4F , - 1) I - ) - I F M + ^ 1 ^ + ^-F 

a a \a) \ 2 a L 2 



a 

\— 

a 



\\2k 1 



(24) 
(25) 



So the equation of state a> DE is given by 

2 (Fj -l)H + 2FjH - (Fj( - 1) (2k/ a 2 ) 



Ude = -1 + 



3(1 - 2Fj)H 2 + (3k I a 2 ) - (1/2) F 



(26) 
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From the data of WMAP + BAO + SN b , the constraints on to DE are io DEQ = -0.999+°;°^ [16]. We 
now define a deceleration parameter 



q = -tt • (2V) 
a 2 



From equations (I22l)and(l23l). we get 



4o = — - 2 = -&DE0 + -z- * -0.592, (28) 

Pmo + Pdeo - 3k/(SnGa 2 ) 2 

here we have used "Q D£0 = 0.728, Q m0 = 0.272" OS from the data of WMAP + 5A<9 + # - The 
critical density of the universe if defined as p c = 3H 2 /SnG. From (1221) it is easy to get 



3// 2 3fc 



PcO = = PmO + PD£0 - - ~ 2 - (29) 



From above, The curvature density is defined as p^ = -3k/(SnGa ). Now we could define 

p c a 2 // 2 

The data of WMAP + 5A<9 + S give the constraints on curvature density, Q k0 = -0.0057^°;°°^ 
[16]. Obviously, It is probable that k > and the universe is close. If the universe is accelerating 
as power-law that a(t) oc t y today, it is easy to get q = y(l - y)/y 2 = -0.592 and y = 2.451. 
We consider the simple case that 

-2A(I,R) = g(I)+f(R). (31) 

From (|2"4"1) and (|25l ), it is easy to get 

SnGp D E = gjl - l -g{I) - l -f(R), (32) 

2k 1 1 

SttGpde = 2 ( gj + Igjjl) H - gJ I - f, R — + -g(I) + -f(R). (33) 

a L 2 2 



The equation of state is 



2{g J + 2g J1 I)H-U2kla 2 ) 

code = -1 H . (34) 

gjI-(l/2)g(I)-a/2)f(R) 



Now we investigate a concrete model -2A(I,R) = aeP R , here R = 6k/a 2 and we assume k > 0. 
The state equation of this model is 

o) m = -\ + 1/3R. (35) 



It is easy to get -0.083 < BR < 0.087 from co DEQ = -0.999™;^ 1160 . The coefficient a can be 
fixed from the equation SnGp DE o = (-l/2)ae^ R °. From (I2TT) and (|22l) . the evolution equation of a 



is 



(I;) 2 = Qk0 (?) 2 + Q ^ o{(ao/fl)2 - 1} + flU (^) 3 + (1 - Omo - Oo - 0*0) (^f . 
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where "Qdeo = 0.73, Q m o = 0.27, Q r o = 0.014 " Rial . The details about the universe evolution 
are shown in Fig. 1 and 2. Obviously in Fig. 1, the evolutions of deceleration parameter q for 
J3Rq = -0.01 or 0.01 are nearly identical to the case /3Rq = which corresponds to the ACDM 
model. Fig. 2 shows that the state function co DE flow to -1 quickly for the cases f3R = -0.01 or 
J3Rq = 0.01 when co DE identically equal to -1 for J3R = 0. 
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Figure 1: The deceleration parameter q - -da /a 2 as a function of a. 
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Figure 2: The parameter code as a function of a. 

Another specific model is -2A(I,R) = ye nl . The corresponding state equation is 

2t]I I+277/ H 



C0 DE = -1 + 



3 \-2qIH 2 



(36) 



When codeo = - 1> H na s two solutions q = and ql = -1/2. When q = 0, it just is a ACDM 
model. The attention will be put on the second case ql = -1/2. The coefficient 7 could be fixed 



6 



from 



(37) 



and y = -3e 1/2 Q. DE0 H 2 . The evolution equation of a is 



Similar to the discussion above, Fig. 3 is the evolution of deceleration parameter q which is similar 
to ACDM model. Fig. 4 shows that the state function a> DE flow to -1 quickly. 




Figure 3: The deceleration parameter q = -aa/a 2 as a function of a. 




Figure 4: The parameter o>de as a function of a. 



The two models discussed above could be regard as the modification to the ACDM model. 
To address the dark energy issues, a changeless cosmological constant is not essential. If the 
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cosmological constant changes slowly and gradually flows to a constant, the observation data 
about dark energy could also be fitted well by this kind of theory. 
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